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Abstract 

We apply the influence-functional method of Feynman and Vernon to the 
study of a single-mode optical field that interacts with an environment at 
zero temperature. Using the coherent-state formalism of the path integral, 
we derive a generalized master equation for the single-mode optical field. 
Our analysis explicitly shows how non-Markovian effects manifest in the ex- 
act decoherence dynamics for different environmental correlation time scales. 
Remarkably, when these are equal to or greater than the time scale for sig- 
nificant change in the system, the interplay between the backaction-induced 
coherent oscillation and the dissipative effect of the environment causes the 
non-Markovian effect to have a significant impact not only on the short-time 
behavior but also on the long-time steady-state behavior of the system. 

Key words: Non-Markovian decoherence, Master equation, Influence 
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1. Introduction 

In realistic conditions, it is impossible to completely isolate a quantum 
system S from its environment E. A proper analysis of the quantum dy- 
namics of S must therefore take into account the decoherence effect of S 
induced by E. Conventional approaches to a theoretical study of the dy- 
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namics of such an open quantum system have treated the interactions be- 
tween S and E perturbatively. Invoking the Born-Markov approximation, 
these yield approximate equations of motion such as the Redfield or mas- 
ter equations {]], 0, 0, 0] for the reduced system we are interested in. The 
Born-Markovian approximation works well when the environment correlation 
;ime te is small compared to the time scale To for significant change in S 
3], 0]. However, in the light of recent experiments (see, for instance, Refs. 
1) @> 0] ) , it is evident that there are many physically relevant situations where 
the Markovian assumption does not hold, and a non-Markovian treatment 
of the open system dynamics is necessary. Furthermore, there is a general 
interest in the fundamental theory of open quantum systems to extend the 
well-developed concepts and methods for Markovian dynamics to the non- 
Markovian case. The development of a general description to open quantum 
system has thus attracted much attention lately. Some recent works explor- 
ing the non-Markovian decoherence properties of quantum systems include 
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Refs. 

In this paper, we consider the exact non-Markovian decoherence dynam- 
ics of a single-mode optical field system - a basic building block in quantum 
communication and computation 17j. Indeed, many of the quantum informa- 



tion processing protocols [l8|jll9j], especially in the field of continuous- variable 



quantum information processing [20J, involve optical fields. In practice, the 
optical field inevitably interacts with surrounding environment, which always 
results in the decoherence of the optical field. Actually, this decoherence ef- 
fect still exists even when the optical field is transmitted in an optical fiber 
2ll . I22J. Decoherence will undoubtedly have a detrimental influence on the 



performance of the protocols. It is thus very important for one to have a 
complete quantitative understanding of the destructive influences of the en- 
vironment. Many of the current quantum optical experiments are performed 
at low temperatures, under vacuum condition. In this case, the main source 
of decoherence is the vacuum fluctuation. Many of the theoretical studies 
on the decoherence dynamics of an optical field to date rely on the Born- 
Markovian approximation 23|, |24], |25|, |26( . We note that, only very recently, 



some phenomenological models on non-Markovian decoherence dynamics of 
optical fields have been studied [l2j, [l3|, [lj]. Based on perturbation, these 



may not capture all the characteristics of the exact non-Markovian dynam- 
ics. A more satisfactory theory, derived from first principles, that describes 
an open optical field system is thus desirable. To this end, we apply the 
influence-functional method of Feynman and Vernon 27], 128] to the study of 
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the system S of a single-mode optical field that interacts with an environment 
E, consisting of a set of harmonic oscillators at zero temperature. 

The Feynman- Vernon influence functional theory allows one to derive the 
exact non-Markovian dynamics of S. In this method, the density operator 
of the combined system, S and E, is expressed as a double path integral. 
The exact dynamics of the reduced system S is then obtained by integrating 
this path integral over the degrees of freedom of E. The effective action 
that governs the evolution of S thus consists of the free action of S and 
an influence functional. All the environmental effects on S are dynamically 
incorporated in the influence functional, and both the backactions from E 
to S and from S to E can be treated self-consistently. In our analysis, we 
mainly address the issues about how the non-Markovian effect takes action 
and what its influence on the decoherence dynamics of the optical field system 
in different parameter regimes is. As an explicit example, we will study the 
time evolution of the so-called Schrodinger cat state [18|, |l9( and see what an 
exact treatment of its decoherence dynamics would yield. 

Indeed, it had been shown that non-Markovian effect shows its significant 
consequence on the decoherence dynamics just by a transient oscillation in 



short time scale [10|, [12|, [13j, [14|]. Here, we will show, besides this transient 
oscillation in short time scale, the non-Markovian effect can also influence 
the behavior of steady state in long-time scale when te — t or te > t . 

Our paper is organized as follows. In Sec. II, we introduce a model of the 
single-mode optical field in an environment, and the coherent-state repre- 
sentation. In Sec. Ill, we present a detailed derivation of the quantum non- 
Markovian master equation using influence functional theory in the coherent- 
state representation. Our generalized master equation has time-dependent 
frequency shift and decay rate, and it reduces to the general Markovian one 
under certain approximation. Sec. IV is devoted to a numerical study of the 
system decoherence dynamics for different environmental correlation time 
scales. In particular, the decoherence dynamics of the Schrodinger cat state 



18l . Il9| is investigated explicitly. Finally, a summary is given in Sec. V. 



2. A single-mode optical field in a quantized radiation field 

We consider a single-mode optical field S that interacted with an environ- 
ment E. The environment, as usual, is modeled by a collection of harmonic 
oscillators. The total Hamiltonian governing the coupled S and E is given 
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by 

H = H S + H E + H I} (1) 

where 

Hs = hw a^a (2) 
is the Hamiltonian of the free single-mode optical field S, 

H E = Ukb\b k (3) 

k 

is the Hamiltonian of the environment E, and 

Hj = hJ2(9ka\ + gtabl) (4) 
k 

is a model for the system-environment interaction, a and a) are the annihila- 
tion and creation operators of the single-mode optical field 5* with frequency 
uj . Similarly, b k and b\ (k = 1, 2, • • •) are the annihilation and creation oper- 
ators of the k-th mode of the environment with frequency u k . The coupling 
strength between S and the kth mode of E is given by g k . The environment 
is assumed to be at zero temperature initially. By this model we will ad- 
dress the decoherence mechanism of S due to its energy exchange with the 
modes of the environmental vacuum fluctuation. We emphasize the model 
and formulation developed in the following are also applicable to many simi- 
lar systems in quantum o ptic s, for example, an optical field in a leaky cavity 



3| or in an optical fiber [21l . |22| , and the noise effect on a nanomechanical 



oscillator 29 



The Hamiltonian H can be expressed, in the coherent-state representation 
30| . in terms of 



H(a, a, z, z) = h{u aa + ^[u k z k z k + (g k az k + g* k z k a)]}, (5) 

k 

where z denotes (zi, z 2 , • ■ •)■ The coherent states \a) = exp(aal~)|0) and 
\z k ) = exp(z k b\.)\O k ) are the eigenstates of a and b k respectively: a\a) = a\a); 
b k \z k ) = z k \z k ). They are nonorthogonal: {a\a') = exp(aa'); (z k \z' k ) = 
exp(z k z' k ). And, they form overcomplete sets: J dfi(a) \a)(a\ = 1 = J dfi(z k )\z k ) (z^l 
with the measures d[i (a) = cx P(~ aa ) d 2 a and dfj,(z k ) = exp ^~ Zfc2fc ) d 2 z k . In the 
path integral approach, one needs to choose a convenient representation. 
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It turns out that employing the coherent-state representation allows very 
straightforward evaluation of the path integrals. In the following section, we 
will use this method to derive the exact decoherence dynamics of the optical 
field system. 



3. Quantum non-Markovian master equation 

3.1. The influence functional in coherent- state representation 

Suppose pr(t) describes the state of our single-mode optical field system 
S plus the environment E as a whole. This total density matrix obeys the 
Liouville-von Neumann equation ihdprit) jdt = [H,pT(t)], which gives the 
formal solution: 

p T (t) = eM~Ht)Pr(0) exp(^ift). (6) 
In the coherent-state representation, pr(t) can be expressed in terms of 

x (a f , z f ; t\a h zf, 0){a h z i |p r (0)|o;-, zj)(a£, z'f, 0\a' f , z' f ; t), 

(7) 

where the resolution of identity has been used. Since we are only interested 
in the dynamics of S, it suffices to work with the reduced density matrix, 
which is obtained by integrating over the environmental variables. This can 
be expressed in terms of 

p(a f ,a' f ;t) = J dp(z f ){a f ,z f \p T (t) \a' f ,z f ) 
= J dp(a i )dp(ai)j(af,ay,t\ai,a' i ;0)p(a i ,a' i ]0). (8) 

Here, we have assumed that the initial total density matrix factors into a 
system part and an environment part, i.e., Pt(0) = p(0) <g> pe(0). Now, 
it remains to determine the effective propagating function for the reduced 
density matrix, 

J{pLf, a'j] t\ai, a-; 0) = J dp(zf)dp(z i )dp(z' j ) 

x (a f ,z f ; t\a h zf, 0}p E (zi, z'f, Q)(a[, z'f, ti\a' p z f ; t). (9) 
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Equation (jH]) contains the forward and backward propagators of the total 
system. These can be expressed as path integrals. To evaluate the forward 
propagator (a/, z/; t\ai, z$; 0), one divides the time interval tf—U into N equal 
subintervals. This is followed by inserting N — 1 copies of the resolution of 
identity, each between a subinterval, and taking the limit of N large. The 
path integral representation of the forward propagator can then be obtained: 

(ckf, Zf; t\ati, Zi, 0} = J D 2 zD 2 aexp(—S[z,z,a,a\), (10) 

with 

S[z, z, a, a] = Ss[a, a] + Si[z, z, a, a] + Se[%, z], (11) 

where Ss, Se, and Si are the (complex) actions corresponding to H$, He, 
and Hj respectively. All the functional integrations are evaluated over paths 
z(r), z(r), a(r), and ot{r) with endpoints z(t) = Zf, z(0) = Zj, a(t) = ctf, 
and a(0) = «j. The backward propagator (c^, z^; 0|«j, zj; t) can be evaluated 
in the same fashion. Substituting Eq. (1101) and a similar expression for the 
backward propagator into Eq. Q we obtain 

J'(af,a'j;t\a i ,a' i ;0) = J D 2 aD 2 a' exp{ — (^[a, a] 

— Sg[a', a'])}J-[a, a, a', a'}, (12) 

where 

J r [a,a,a',a'] = J dfi(zf)dfi(z i )dfi(z' i )D 2 zD 2 z 

i 

x p E {zi, z-; 0) exp{-(S E [z, z] + Si[z, z, a, a] 

-^[z',z']-5;[z',z',a',a'])}, (13) 

is the influence functional containing all the environmental effects on S. 

3.2. Evaluation of the influence functional and effective propagating function 

Now, we calculate explicitly the influence functional of our model. Using 
the Feynman's procedure, one can obtain the path integral representation of 
the forward propagator, Eq. ljTOl) . for our system with the component actions 
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as 



Ss[ct, ot] = —ihaa (t) + / dr[ihaa(T) — Hs(a, «)]}, 

Jo 

S E [z, z] = y^{-ihz k z k (t) + / dr[ihz k z k (T) - H E (z, z)]}, 
k Jo 

Sj[z, z, a, a] — — / drHi(a, a, z, z). (14) 
Jo 

The path integral with respect to the environmental variables z can be evalu- 
ated by the saddle point method under the boundary conditions -2^(0) = z k i, 
z k {t) = z k f. We have the equations of motion as 

z k + iu k z k = -ig* k a, z~ k - iu k z k = ig k a, (15) 

where a and a are treated as external sources. Substituting the solution of 
Eqs. ( 1151) into Eq. ffTU]) . one can determine the desired path integral. It 
is noted that the prefactor under the contribution of stationary path in the 
coherent-state representation is equal to one, and the saddle point approach 
to the evaluation of the environmental part here is exact 31] • The path inte- 
gral with respect to the environmental variables z', in Eq. (fl3l) . can similarly 
be obtained. As explained in our introduction, we take the environment to 
be at zero temperature, i.e. p^(zj,z^;0) = 1. Together with the results of 
the above path integrals, Eq. (1T5I) yields the influence functional 

J r [a,a,a',a'} = exp{J dr J dr'[p,(r — r')(a' 
-a){r)a{r') + /x*(r - r')a'{r'){a - a')(r)}}, (16) 



with /i(x) = ^2 k e luJhX \9k\ ■ Here, we have repeatedly used the Gaussian 
integral identity / ^± e -^+^f^) = ±/(£). 

In the derivation of the influence functional above, we have treated both 
the backactions of the environment on the system and the system on the 
environment self-consistently. All the effects of the environment E on S 
are collected in the influence functional, which results in a correction to the 
action of the free single-mode optical field 5", 

J{pL^ a'j] t\oii, c^; 0) = J D 2 aD 2 a exp{aa(t) 
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+ a'a'(t)— / dr[aa + a'a' + iHs(a,a) 
Jo 

— Hs(at', a')]} J- [a, a, a', a']. (17) 

To evaluate the path integral in Eq. (fT7]) . we again employ the saddle point 
method and obtain the two equations of motion: 



= a + iuj^a + / dr'fi (r — r') a (r) , 
Jo 

= a-iu a'+ [ dr' // (r - r') a' (r') , (18) 
Jo 

with the boundary conditions a (0) = a! (0) = a[. The integro-differential 
equations render the reduced dynamics non-Markovian, with the memory 
effect of the environment registered in the kernel that is nonlocal in time. 
The solution of the integro-differential equations (ITS]) can be expressed in 
terms of a complex function u{t) as 

a(r) = ajti(r), «'(t) = o£u(t), (19) 

with the boundary condition w(0) = 1. Substituting Eqs. (fTDl into Eq. (fiTI) 
and using Eqs. f|T8l) . we obtain the expression of the effective propagating 
function of the reduced system as 

J(oif, a.y, t\d(i, a[; 0) = exp{[ua>fOii + ua^a'^ 

+ (1 - M 2 )c^]}, (20) 

where the dependence of u on time is not shown explicitly for abbreviation. 

3.3. The non-Markovian master equation 

Now we can derive the master equation by computing the time derivative 
of Eq. ([8]). First, from Eq. (120]) . we can write down the following identities 



1SJ _, 1SJ 
oiiJ = -t— ; ®iJ = -77, (21) 



which will be used to remove from the time derivative of J its dependence 
on a.i and a\. After taking time derivative to Eq. (jSJ) and substituting Eqs. 
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Figure 1: Comparison of the decay rate T(t) (left) and the frequency shift 5u(t) (middle), 
and the purity (right) of the Schrodinger cat state between the non-Markovian (solid 
line) and Markovian (dashed line) results in the weak coupling and short environmental 
correlation time regime. The parameters uj c /ui n — 50.0, 77 = 0.1, and f3 = 1.0 are used in 
the numerical calculation. 




Figure 2: Comparison of the decay rate T(t) (left), frequency shift 6u>(t) (middle), and 
purity p(t) (right) of the Schrodinger cat state between the non-Markovian (solid line) and 
Markovian (dashed line) results in the strong coupling and long environmental correlation 
time regime. The parameters uj c /ujo — 1.0, f] — 5.0, and (3q = 1.0 are used in the numerical 
calculation. 




Figure 3: Comparison of the decay rate T(t) (left), frequency shift 5u(t) (middle), and 
purity p{t) (right) of the Schrodinger cat state between the non-Markovian (solid line) and 
Markovian (dashed line) results when te ^> t . The parameters oj c /^o = 0.2, 77 = 5.0, 
and /3q = 1.0 are used in the numerical calculation. 
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(I2"0~j) and (T2~Tj) into it, we obtain the evolution equation 

■(- ' +\ r -rv+u-^"'"' 5 *) Sp(6c,a';t) 

p(a, a;t) = {— lillt) [a — ~ a 

oa oa 

+r(t)[2 fcfe ~ a ^a 6a~ a]h (22) 



where 



- = -T(t)-in(t), (23) 



and from here on we drop the subscript /. 

Next, we introduce the following functional differential relations in the 
coherent-state representation: 

_Sp(a,a';t) t 5p(a,a';t) t 

a. — < > a'ap{t), a < > p{t)a'a, 

oa oa 

5 2 p(d, a':t) , . + /n , s 

" (24) 

oaoa 

with which we arrive at our final operator form of the non-Markovian master 
equation 

p{t) = ~[H\t),p(f)] +Y{t)\2ap{t)a ] - a ] ap{t) - p{t)a ] a] } (25) 

where H'(t) = hVl{t)a)a. This is the exact master equation for the reduced 
system dynamics. Our exact master equation (1251) is similar to the non- 
Markovian master equation of a two-level atom in a bosonic environment 



32], |4|. This similarity for the fully different systems comes from the fact 
that both the derivations are based on the rotating wave approximation and 
the initial vacuum state of the environment. O(t), named as time-dependent 
Lamb shifted frequency in the two-level atom system, plays here the role 
of a time-dependent shifted frequency of the single-mode optical field S in- 
duced by the environment E. T(t) represents a time-dependent decay rate 
of the optical field. We emphasize that our derivation of the master equa- 
tion goes beyond the Born-Markovian approximation and contains all the 
backactions between E and S self-consistently. It can be seen that the only 
difference between Eq. (1251) and the master equation under Born-Markovian 
approximation [3| is the time-dependent coefficients. So we argue that all 
the non-Markovian characteristics reside in the time-dependent coefficients 
of the generalized master equation. 
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The time-dependent coefficients in our generalized master equation, de- 
termined by Eq. (I2"5j) . essentially depend on the so-called spectral density, 
which characterizes the coupling strength of the environment to the sys- 
tem with respect to the environmental frequencies. It is defined as J(oo) = 
Ylk \9 k \ 2 ~ Wk )- m ^ ne continuous limit of the environmental frequencies, 
we model in our work the coupling of the optical field with the environment 
has a spectral density as 

J(u)=T}u[ — ) e"~, (26) 
where 77 is a dimensionless coupling constant, and uj c is an ex pon ential cutoff 



frequency. The environment is classified into three categories |33|: sub-Ohmic 



if < n < 1, Ohmic if n = 1, and super-Ohmic if n > 1. Different spectral 



densities manifest different non-Markovian decoherence dynamics [33j. The 
spectral density Eq. ( 1261) is motivated by the physical consideration that not 
all the modes of the environment give the same contribution to the coupling 
to the optical field system. Thus, the spectral density form is physically 
reasonable to model the decoherence of our optical field system. Actually 
such spectral density, especially the Ohmic one, is widely used in the the 
decoherence analysis of optical fields in the scenario of the continuous variable 
quantum information processing 



Before presenting our numerical results in the next section, we show how 
our generalized master equation reduces to the conventional master equa- 
tion by introducing the relevant Markovian approximation. By defining new 
dynamical variables as x(t) = a(r)e' luJoT , we can recast the first equation of 
Eqs.fllH]) into 

POO pT 

±{r) + / dujj{uj) / dr'e^-^-^xir') = 0. (27) 



Then, invoking the Markovian approximation, x(t') ~ x(t), namely, approx- 
imately taking the dynamical variable to be the one that depends only on the 
present time so that any memory effect is ignored. The Markovian approxi- 
mation is mainly based on the physical assumption that the correlation time 
of the environment is much smaller compared with the typical time scale of 
the system evolution. Also under this assumption we can extend the upper 
limit of the r' integration in Eq. (j2"7|) to infinity and use the equality 

lim / d r ' e ±*(-o--)(T-r') = _ ^ TlV f — 1 — \ ^ ( 2g ) 
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where V denotes the Cauchy principal value. The integro-differential equa- 
tion (12T1) thus reduces to a linear ordinary differential equation. The solution 
of x(t), as well as ck(t) can then be obtained readily, which results in 

U (t) = e -i(wo-««)T-7rJ(«o)T ) (29) 

with 5ui = V J °° • Using this solution, one can verify from Eq. fl23l) 

that 

r(t) = 7tJ(cj ), fi(t) = c; - (Jw, (30) 

which are precisely the coefficients in the Markovian master equation of the 
optical system 0]. 

In the next section, for definiteness, we consider Ohmic environment E. 
The characteristic time scale te of the environmental correlation function in 
the Ohmic case is roughly inversely proportional to the cutoff frequency uo c 
in Eq. (I26I) . i.e., te — l/co c [111]. The cutoff frequency cu c , which is originally 
introduced to eliminate infinities in frequency integrations, therefore also 
determines if the dynamics of open system S is Markovian or non-Markovian. 
Our non-perturbatively derived exact results can allow us to explore all these 
possibilities. 



4. Numerical results and discussions 

To illuminate the non-Markovian decoherence dynamics of S, we consider 
the following initial state of the optical field: 

p(0) = ^[|A)XA)I + I - A)X-A)| + IA))(-A)| + I - A>)(A>|], (31) 

where N = 2(e' /3 °' + e - '*' ) is a normalization constant. This is known as 



the Schrodinger cat state and has been produced experimentally [la . Il9 



After some straightforward calculations, we obtain, via Eqs. (|HJ) and (|20|) . 

Pit) = ^[e m2 - m \\m\ + \-P){-m 

+e -(l*| 2 -|/3| 2 ) ( |_^ ) ^| + | /3)( _^| )]; (32) 

where (3 = 0Qu{t). From Eq. (1521) . the purity which is defined as p(t) = 
Trp 2 {t) can be calculated readily as 

p(t) = l_i e 2\M 2 + e -2\M 2 + e 2\M 2 -m 2 

+e -m\ 2 +m 2 +A y (33) 
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In Fig. [T], we plot the numerical results of the decay rate T(t), frequency 
shift 6u)(t), and purity pit) when t e <C r . We note that the exact T(t) and 
8u(t) differ from their corresponding values obtained via a Markovian treat- 
ment only over a very short time interval. For t < te, both coefficients grow 
very quickly, while for t > te, they gradually approach the corresponding 
Markovian results, given by Eqs. fl3"Ul) . as t approaches r . The finite, almost 
steady, positive decay rate guarantees the irreversibility of the system dy- 
namics. Accordingly, the exact time evolution of the purity shows only slight 
deviation from the Markovian results. The Schrodinger cat state eventually 
evolves to a steady state, namely the ground state of S: p g = |0)(0|. Clearly, 
in this case, the backaction of the environment E has negligible effect on the 
dynamics of system S, and the Markovian approximation is applicable. We 
say the system dynamics is mainly governed by the dissipative effect of E. 

Figure [2] shows the decay rate T(t), frequency shift Su(t), and purity 
p(t), when te = To- In this case, the backaction of the environment E has 
a considerable impact on the dynamics of our system S, and the Marko- 
vian approximation is not applicable. Firstly, we note that, in contrast to 
the Markovian treatment, the decay rate can take negative values. Physi- 
cally, this corresponds to S reabsorbing a photon from E, which will lead 
to an increase in the photon number of S [4|. Next, more interestingly, we 
note that the decay rate approaches zero asymptotically, which dramatically 
differs from the Markovian results. Consequently, the exact evolution of 
the Schrodinger cat state is drastically different from the Markovian results. 
In particular, it eventually evolves to some steady state, which is not the 
ground state p g . The backaction of the environment causes the system to 
undergo transient oscillations, which is characteristic of non-Markovian dy- 
namics. From previous studies 12|, [l3|, ll4| , one would have concluded that 
non-Markovian effects only show up in short-time dynamics. Our results, 
however, clearly show on the contrary that non-Markovian effects can also 
have an influence on the long-time behavior of the system dynamics and the 
final steady state of S. This is because that the dissipative effect of E is 
strongly counteracted by the effect due to the backaction of E. 

The results under the condition t e ^> r is shown in Fig. [3J Because of 
the extremely long memory effect of the environment, the backaction from 
E to S is so strong that it governs the dynamics of S. This causes the decay 
rate to oscillate over a very long duration. This is typical of non-Markovian 
dynamics The 'equilibrium' position for the oscillation of T(t) is not at 



zero, but a small positive value. The positivity of the equilibrium position 
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of the decay rate oscillation makes the system dynamics experiences weak 
dissipation. Such weak dissipation is verified by the time evolution of the 
purity of the Schrddinger cat state in Fig. [31 The evolution of p(t) also 
shows that the coherent oscillation induced by the backaction of E persists 
for a very long time, even as the state approaches the ground state. 

In summary, depending on Te in comparison to r , the decoherence dy- 
namics of S shows different behaviors. For < t 0i the exact results show 
only slight deviation from the Markovian ones. Since the effect due to back- 
action of E is extremely small, the system dynamics is mainly governed by 
the dissipative effect of E, and the widely used Markovian approximation is 
applicable. For r E = r , the considerable backaction induced by the near 
resonant interaction between E and S counteracts the dissipative effect, and 
results in zero decay rate asymptotically. This causes dissipation to cease 
before the system decays to its ground state. That is, the steady state is 
no longer the ground state, unlike the previous case. When te ^> t , the 
backaction of E is so strong that it governs the dynamics of S. The decay 
rate of the system oscillates about some equilibrium position over a very long 
period of time. The positivity of such an equilibrium position guarantees the 
overall weak dissipation effect of E. 

5. Conclusion 

In this work, using the influence-functional method of Feynman and Ver- 
non, we investigate the exact decoherence dynamics of a single-mode optical 
field S in an environment at zero temperature. We derive an exact gen- 
eralized master equation for S. The equation's time-dependent coefficients 
depend on the environmental spectral density. These determine the exact 
dynamics of S. We conclude from our numerical results that E exerts two 
competing influences on S. One effect, Si, is dissipation, and is responsible 
for the decoherence of S. The other, £ 2 , is due to the backactions of E on 
S. The backaction of E on S means that E with its state changed due to 
interaction with S, in turn, exerts its dynamical influence back on S. This is 
the physical origin of the non-Markovian dynamics of S. In the conventional 
Born-Markovian treatment, one generally neglects the backaction of the en- 
vironment. So, in that case, we do not see the effect due to backaction on the 
dynamics of S. The degree of manifestation of £\ and £2 in the dynamics of 
S depends on te in comparison with r . For te <C t , £\ dominates and £2 
only gives rise to a transient coherent oscillation of S. The state of S evolves 
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to the ground state p g , which is coincident with the Markovian result. If 
t e — T o, £1 and E-2 are comparable and their effects counteract each other. 
The state of S asymptotically evolves to some steady state, which is not the 
ground state of S. Finally when te ^> tq, £2 dominates and governs the 
dynamics of S. The state of S eventually approaches to the ground state but 
never quite reach it. 

The theory we have established is a non-perturbative description of the 
exact decoherence dynamics of a single-mode optical field. It is a general- 
ization of the well-developed Born-Markovian treatment of such system. It 
can serve as a useful basic theoretical model in analyzing the non-Markovian 
decoherence dynamics of optical fields employed in practical quantum in- 
formation processing schemes. It should be noted that although only the 
Ohmic spectral density is considered here, it is straightforward to generalize 
our discussion to the non-Ohmic cases. 
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